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Discovery of a high-redshift Einstein ring* 
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Abstract. We report the discovery of a partial Einstein ring of radius 1'.'48 produced by a massive (and seemingly isolated) 
elliptical galaxy. The spectroscopic follow-up at the VLT reveals a 2L t galaxy at z = 0.986, which is lensing a post-starburst 
galaxy at z = 3.773. This unique configuration yields a very precise measure of the mass of the lens within the Einstein radius, 
(8.3 ± 0.4) 10 1 1 /i^q Mq. The fundamental plane relation indicates an evolution rate of dlog(M/L)g/dz = -0.57 ± 0.04, similar 
to other massive ellipticals at this redshift. The source galaxy shows strong interstellar absorption lines indicative of large gas- 
phase metallicities, with fading stellar populations after a burst. Higher resolution spectra and imaging will allow the detailed 
study of an unbiased representative of the galaxy population when the universe was just 12% of its current age. 
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1. Introduction 

One of the key issues in galaxy formation within the current 
ACDM framework of structure formation is the mass assem- 
bly histories of galactic halos. In this paradigm, the growth of 
halo mass through mergers produces giant galaxies and star for- 
mation appears rather late during this process. Measuring the 
evolution of the mass-to-light ratio hence constrains directly 
this formation scenario, and provides clues on the evolution of 
the fundamental plane. Various deep surveys have uncovered 
different galaxy populations, but the selection cr iteria produce 
biased sample s : UV-selected JSteidel et al.l2003l) and narrow- 
band selected Jhu" et all2002l) samples are sensitive to actively 
star-forming galaxies and biase d against quiesce nt, evolved 
systems, while sub-millimeter and near- 

infrared surveys ( McCarthy et al .1 1200 it lAbraham et al.ll2004T) 
select dusty starburst galaxies and very red galaxies respec- 
tively. It is not clear how representative these samples are of 
the population of distant galaxies as a whole and how they are 
related to present-day galaxies. In contrast, selection through 
gravitational lensing is not biased, as any distant galaxy which 
happens to have a massive deflector along the line of sight 
can be amplified. This technique has yielded many examples 
of very distant galaxies, amplified by the mass distribution of 



foregroun d clusters of g alaxies. ( Ellis et alJl200lt IPettini et alJ 
2002| IPello et afl Eo04). Configurations where lenses are iso- 
lated and distant (z ~ 1) galaxies are much rarer and yet pro- 
vide direct measures of the total mass without any assump- 
tions on stellar evolution. Whilst the total mass enclosed within 
the projected Einstein radius 0e can be well measured with 
systems presenting a few (~2-4) images of the lensed galaxy, 
many more constraints a re given by the observa tions of nearly- 
complete Einstein rings jKochanek et al.l200lh . as they do not 
suffer from the well-known ellipticity-shear degeneracy due to 
the many data along different position angles. 

Although many arcs have been discovered, associated with 
massive galaxy clusters and their domi nant central galax- 
ies, th ere are very few optical rings ( Miralda-Escud e~& LehaJ 
119921) or arcs, especially at high redshifts. Only a few sys- 
tems with lenses above zi ~ 0. 9 are known : CFRS 
03.1077 fa = 0.938,z s = 2.941. ICrampton et alJ | 2002h 
MG 0414+0534 fa _ .958, z s = 2.62. iHewitt et alJll992l 
iTonrv & Kochanekll 19991). MG 2016+112 fa = 1.004 z s = 
3.273. Lawrence et all 1 19841: iKoopmans & Treul tZ002b. and 
possib ly J100424.9+122922 fa ~ 0.95, z s = 2.65.lLacv etal 
120021) and GDS J033206-274729 fa ~ 0.96. lFassnacht etal 



* Based on observations carried out at the ESO/VLT under pro- 
grams 68A-01706(A), 71A-0102(A) and 72.A-01 11(A) and at La 
Silla 3.6-m telescope under DDT. 



2004). In all cases the arcs cover less than 60 degrees around 
the central lens and hence are not nearly-comple te Einstei n 
rings. The only optical ring reported so far ( Warren et alJl9 96). 
over some 170° with a radius 1'.'08, is a galaxy at z s = 3.595 
lensed by an elliptical galaxy at a much smaller zi = 0.485. 
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Fig. 1. (a) Image of the ring with 075 seeing, (b) best configu- 
ration for the lens and the source, with the position of the caus- 
tics (cf. text for parameter definitions), (c) best-fitting image 
convolved with the seeing, and (d) the predicted image at the 
resolution of HST/ACS. 



the fit was found in the 10-dimensional parameter space, we 
coupled a state-of-the-art genetic algorithm developed by one 
of us with the gravlens code (Kee t on! 1200 II) . Briefly stated, 
the parameters are coded in the chromosome of an individual 
which can produce two offsprings from its coupling with an- 
other individual (Charbonneau 1995). For each individual in 
each generation (i.e. one configuration in parameter space), we 
compute the resulting image, convolve it with the seeing of the 
VLT image, and measure its likelihood. These steps are illus- 
trated in Fig. ^ After a few thousand generations the genetic 
algorithm finds the absolute maximum of the likelihood, and 
several tests were made to assess that different runs always con- 
verged to the same solution. From the ensemble of generations 
and individuals the likelihood contours can be computed and 
the partial correlations and degeneracies between parameters 
explored. The procedure will be presented in detail elsewhere. 



Following iKassiola & Kovnerl dl993h : iKormann etail dl994 : 
iKeeton et alJ ( 1 19981) . a projected surface mass density for an 
isothermal ellipsoid (SIE) written in polar coordinates (r, <p) 
and expressed in units of the critical density for lensing, 



K(r,4>) = 



b_ 
Yr 



1 + q 2 



(1 +q 2 )-(l -q 2 )cos24>) 



Here we report the discovery of a fourth confirmed system, 
dubbed FOR J0332-3557 (03 1 ?32™59 s :94, -35°57'5 177, J2000), 
in a sight line through the outskirts of the Forna x cluster, w here 
the re ddening of our Galaxy is E(B - V) = 0.02 JSchlegel et alJ 
Il998l) . This new system is remarkable on two accounts. First, 
it is a bright, almost complete Einstein ring of radius l'/48, 
covering some 260° around the lens, extending over 4"x3" 
(§2) and with a total apparent magnitude R c = 22.2. As dis- 
cussed in §3, the lens appears to be an elliptical galaxy at 
Z\ = 0.986. Second, the lensed source is a galaxy at redshift 
z s = 3.773 (§4). A flat FRW metric with Q. m = 0.3, f2 A = 0.7 
and Ho = 70/i7o km s _1 Mpc~'is assumed. 



2. Modeling the gravitational lensing system 

Full details of the observational photometric and spectroscopic 
setups will be reported elsewhere (Cabanac et al. 2005). Given 
the excellent image quality in the R c band, with 075 seeing, a 
detailed modeling of the lensing syste m can be carried out. We 
followed the procedure pioneered by Kochan ek et alJ (|2001) 
with two important differences. First, not only we use the de- 
tailed shape of the ring (ridge) and the flux along the ridge, but 
also the width of the ring along each radial direction. These 
three sets of observables were incorporated into a likelihood 
function. Second, since we do not want to impose any a priori 
information, the problem has a high dimensionality: 5 param- 
eters to describe the properties of the lens (given a mass pro- 
file: the normalisation, ellipticity and orientation, plus an ex- 
ternal shear and its orientation) and 5 for the source (position 
with respect to the lens, effective radius, ellipticity and orien- 
tation). Hence, to guarantee that the maximum likelihood of 



where q < 1 is the axial ratio. The radius r and the parameter b 
both have dimensions of length (expressed here in arcsec). b is 
equ al to the E instein radius 6e in the spherical case, Following 
iHutereret al] 12004) we use an Einstein radius 

6» E = Z?exp[(0.89e) 3 ], 

for a non-spherical SIE of ellipticity e — 1 - q. 

Figure [2 shows the best-fit example with b = 17318 + 0.02 
and e = 0.55 + 0.01, yielding 6» E = 1748 + 0702, that is, 
11.8 kpc, with possibly a systematic error of 9% due to 
the different possibl e normalisations of isothermal ellipsoids 
dHuterer et al)l2004l) . Given the redshifts of the lens (§3) and 
the source (§4), the total mass within the Einstein radius is 
M E = (8.3 + O.2)lO ll /i7 o 1 M . The lens is 20% less mas- 
sive than the elliptical le ns in MG 2016+112 at z\ = 1.004 
jKoopmans & TrevJ2002l) . Quite independently of the assumed 
elliptical potential, the derived magnification of the source is 
12.9. The mass ellipticity appears to be larger than the light 
ellip ticity of 0.2 (§ 3) but this is often found in other systems 
jKeeton et alll998t) . 

The large number of constraints provided by the ring lifts the 
degeneracy between the flattening of the potential and the ex- 
ternal shear. For a wide range of potentials we consistently find 
high flattenings and small shears. The small value of the shear 
(Fig. lb) y = 0.08 + 0.01 is typical of those found in loose 
groups of galaxies. Given its direction (48°) it is possible that 
the shear is caused by a faint nearby galaxy which lies £ = 577 
away at a position angle 38° ± 0°5 and with a photometric red- 
shift of ~1.27. If this is the case, the galaxy would have a SIS 
Einstein radius of 8e= 2y£ = 0791 and a projected velocity dis- 
persion of 260 km s 1 . It would imply a massive galaxy, incon- 
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Table 1. Photometric properties of the lens galaxy 



Band 


effective radius R c 




>' 


Mtot 




["] 


[^70 kpc] 


[mag-arcsec 2 ] 


[mag] 


R c 


0.36+0.05 


2.9+0.4 


22.20= 


tO.10 


22.42 


i 


0.36 


2.9 


20.70= 


t0.20 


20.92 


.1 


0.41+0.01 


3.3+0.1 


18.64= 


h0.03 


18.58 


H 


0.43±0.01 


3.4±0.1 


17.70= 


h0.03 


17.54 


K s 


0.50+0.01 


4.0+0.1 


17.13= 


t0.03 


16.63 



1 < /j > is the average surface brightness, not to be confused with the 
surface brightness at the effective radius yu(/? e ) = yu c . 

sistent with its observed brightness. On the other hand, the R c 
and B images show a diffuse population of blue galaxies within 
one arcmin of the Einstein ring, and the near-IR images in J, 
H, and K s provide evidence for a group of galaxies possibly 
associated with the lens or else with a proto- cluster at red shift 

2 as the photometric redshifts with HYPERZ (iBolzonella et alJ 
2000) show a bimodal probability distribution function with 
peaks at both redshifts. There is no evidence, however, for a 
red sequence at z ~ 1 . 



3. The lens galaxy 

The BVRJJHK S multicolour images obtained at ESO/VLT/La 
Silla and the optical spectrum of the lens (Figure [3 constrain 
the spectral energy distribution (SED) unambiguously and al- 
low the precise identification of the lens as a quiescent galaxy 
at redshift zi = 0.9 86 + 0.005 

We used GALFIT JPeng et al.ll2002l) to fit the profile of the 
lens, using a PSF-corrected profile with De Vaucouleur, Sersic 
and exponential parameterizations, masking the ring and us- 
ing a range of sky level (average±3cr). This yields best-fit de 
Vaucouleur 1 profiles with effective radii and average surface 
brightnesses given Table 1, and consistently small ellipticities 
e — 1 - q — 0.20 ±0.1 in all bands. This distant lens appears 
rounder than the nearb y ones, where the mea n observed light 
ellipticity is 0.31 ±0.1 8 Jjoreensen et all 1995b . 
The fitted SED (Fig. El allow s one t o infer properly the k cor- 
rections following PoggiantiJ dl997l) : we find B rest - R co \, s = 
-0.65, B rest - iobs = 0.73 and Z? rest - John = 3.08, yielding a con- 
sistent Z? rest = 21.7 ± 0.1. Using other templates which provide 
similar good fits to the SED yield the same result. For sake of 
completeness, it should be noted that HYPERZ can also fit (with 
a lower probability) bluer templates of you nger age s with ex- 
tinctions of up to A v = 0.8 mag, assuming a Seaton extinction 
law (1979). The definitive internal extinction of the lens can 
only be better constrained with a higher quality spectrum of 
higher resolution. 

If we assume that the lens is an elliptical galaxy, we can com- 
pare its measured optical properties with the nearby and high- 
redshift sample. The rest-frame Mb = -22.3 is twice brighter 

1 Sersic profiles also fit the surface brightness with values of n > 3, 
the dominant parameter in the Chi-square fit is the background sub- 
traction, but since we find no hint of a disk in the R band (blue rest- 
frame), we used de Vaucouleur profile parameterization 
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Fig. 2. Top : The spectrum of the lens (middle black line), the 
underlying sky emission (bottom red line), and the HYPERZ - 
derived best-fit synthetic spectrum (top blue line). Bottom : 
Observed BVRJJHK S photometry and the spectral energy dis- 
tribution of the best-fitting spectrum shown above. 



than the typical galaxies at z ~ 1 which have Mg = -21 .5 ±0.2, 
and its rest-frame U—V colour is redder (1.31 ) than the averag e 
L* galaxy at this redshift (U- V = 1.10 ±0.1 . Bell et alJ2004 . 
while its rest- fr ame U - B = 0.18 a ppears bluer than similar 
field ellipticals ddebhardt et al.l2003l) . 

The rest-frame absolute B-band surface brightness SBb, in- 
ferred from the R c (resp. J) images, yield SBb = 18.57 ± 0.10 
(18.75) mag arcsec~ 2 after the (1 + z) 4 cosmological correc- 
tions. The fundamental plane relation of elliptical galaxies can 
be written as 

n = log(fl e /kpc) - a logo- -B SB B , 

where R e is the effective radius, <x is the velocity dispersion, 
and we adopt a = 1. 25, _ 0.32 and To - -9.04 - log 
hjo (Treu & Koopmans lTreu & Koopmanslll2003l) '). Taking the 
velocity dispersion inferred from the lens modeling (§2) as 
cr = 302 ± 12(syst) km s _1 , the evolution in the mass to light 
ratio becomes 

Alog(M/L B ) = -(Tj - r )/2.5j8 = -0.57 ± 0.04 

where the uncertainty takes into account the systematics in 
the dispersion velocity and the effective radius. This evo- 
lution rate is cons istent with what is observed in field or 
cluster ellipticals JGebhardt et al J l2003t iRusin et alJ [2003; 
iKoopmans & Treu l|2003l) , and w ithin the scatter seen for mas- 
sive galaxies jTreuetaDl2005h at these redshifts. It also im- 
plies a relatively large formation redshift, in a simplified sce- 
nario where star formation proceeds in a single burst. This is 
consistent with the SED-inferred stellar population. We note 
that the radio ring MG 11 31 +0456 (whose lens is at z\ =0.844, 
iTonry & Kochaneldb oOO) also yields a smaller evolutionary 
rate jRusin et al.l2003l) . 
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Fig. 3. Top : The spectrum of the source galaxy, showing 
strong lines from the interstellar medium, yields a redshift of 
3.773±0.003. Bottom: The shape of the continuum and the 
depth of the lines are consistent with a population of stars 
having past a major episode of star formation, very similar to 
the post-starburst galaxies shown in comparison, both distant 
(cB58) or nearby (NGC 1741). The sky emission and some tel- 
luric lines are also indicated. 

Assuming that the evolution of the stellar mass-to-light ratio 
is the same as the effective one derived above, and taking a 
value of 7.8 + 2.2 IiioMq/Lbo for the local galaxies, results 
in (M»/Lb) = 2.1 + 0.1 hjoM Q /LBe, entirely consistent with 
relatively old stellar populations (as corroborated by the best 
fitting SED, Fig. [2J. However, at Einstein radius 6 E the ratio 
M(< E )/L B (< 0b) = 7.1 + 0.5 h 70 MJL BQ . 



4. The source galaxy 

The VLT/FORS low-resolution spectrum of the lensed source 
(Figure^} reveals a galaxy at a redshift z s = 3.773 + 0.003. Its 
apparent total magnitude R c = 22.6 makes it about one magni- 
tude brighter tha n the brightest Lym an-break galaxies (LBGs) 
found at z ~ 3.5 ( Steidel et al. l2003h . The source is sufficiently 
bright to make future follow-up kinematic and abundance stud- 
ies of its interstellar medium sim ilar to what was p ossible on 
the z s = 2.73 lensed galaxy cB58 JPettini et all2002l) . 

The overall spectrum comp ares well with the h igh luminos- 
ity Lyman-break galaxies ( Shapl ev et al.l 2003) which have 
Lyq in absorption at z ~ 3 or even with SDSS Jl 147-0250 
( Ben tz et alJl2004l) at a lower redshift. The continuum appears 
to be rather flat, bluer than cB58 but redder than NGC 1741, 
and it is remarkable that there are no emission lines (with the 
possible exception of C m] AA 190.7,190.9). It is clearly not 
like the H n g alaxy detected at a s imilar redshift (z = 3.357) in 
the Lynx arc jFosburv et alJl2003i) . The flat continuum and the 
absence of Lyq in emission, as seen in many nearby starbursts 



and about half the LBGs, is consistent with a post-burst stel- 
lar population, where the absorption stella r Lyq line be comes 
important after a few million years ( Valls-Gabaud 1993) inde- 
pendently of metallicity or extinction. 

Although the resolution of the discovery spectrum is not high 
enough to make a kinematic study of the different lines, the 
signal-to-noise is sufficient to allow the identification of many 
photospheric and interstellar medium absorption features with 
the following rest-frame equivalent widths (in nm) : C n 
A133.5 (0.15+0.02), C iv ,L1 154.8, 155.1 (0.08+0.02), O i/Si n 
[blend],Ll 130.2,130.4 (0.17±0.02), Si n ,1126.0 (0.30+0.04), 
Si ii ,1152.7 (0.25±0.04), Si iv ,1139.4 (0. 13+0.02), Si iv ,1140.3 
(0.08+0.02). There are also tentative detections of other impor- 
tant lines, such as Fe n ,1160.8, C m ,1117.6, and Al n ,1167.0, 
along with possible features produced by P v and O i. The 
metallicity index at 143. 5nm is contaminated by sky lines, but 
the pattern that emerges from these values, although prelimi- 
nary, is that the gas-phase abundances of these interstellar ab- 
sorpti on lines are r ather high, especially in comparison with 
cB58 JPettini et alfclioll) . In this respect, the source galaxy ap- 
pears more similar to the LBGs at z ~ 5 JAndo et al] 12004) 
where the ISM lines are stronger than at z = 3 (with the no- 
table exception of carbon). These indications point to a very 
rapid enrichment by type II supernovae, associated to bursts of 
star formation, in many respects similar to the patte rn seen in 
star-fo rming galaxies at z ~ 2.2 (Shapley et al IShaplev et alJ 
il2004 v 

The flux at 150 nm has been traditionally used to measure the 
star formation rate from UV spectra. Calibrating on the well- 
measured R c magnitude, we derive a flux density at 150 nm of 
2.15 fiJy (assuming a conversion factor of 2875 Jy for a R c = 
source), which translates into SFR nv » 31 (A/12.9 ) -1 M Q 
a -1 , adopting the standard relation (K ennicutl 19981 and where 
A is the gravitational amplification produced by the lens (§2). 
There are two important caveats. First, this could be a lower 
value because we have neglected the dust absorption, both in- 
ternal to the source galaxy and in the lens galaxy. Second, the 
relation only holds for a constant star formation rate, while it 
seems more likely that star formation proceeded in a series of 
bursts. Catching the galaxy after the burst has finished, or is 
decreasing, implies that the UV luminosity is a strong func- 
tion of age and cannot be properly related to the rate that pro- 
duced these stars unless their age can be measured. All in all, 
the large inferred rate is better used for comparison with other 
galaxies, in which case this source appears to have been ex- 
tremely active. Although it is more than an order of magnitude 
below the extremely bright SDSS Jl 147-0250 JBentz et alJ 
2004) it is entirely con sistent with the rates inferred for LBGs 
at 7 ~ 3 JShaplev et all2003l) or the ~ 40 M a" 1 rate of cB58 
dPettini etal.l2 002). and would explain the rapid metal enrich- 
ment of its interstellar medium sketched above. 

The modeling of the lensing system (§1) also yields con- 
straints on the effective size of the source galaxy, 6* e ff = 0'.'304 
(2.16 kpc), assuming an elliptical gaussian source 2 . Other 

2 the actual shape and position angle of the source are not strongly 
constrained by present images 
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elliptical shapes yield similar values. The compactness of the 
region which gives rise to this emission, combined with its ten- 
tative high metallicity may perhaps be associated with a pro- 
genitor of a present-day bulge. 

Additional observations are clearly required to further con- 
strain this lens, and to study the interstellar medium and the 
stellar populations at a look-back time of 88% of the present 
age of the universe. 

Acknowledgements. We are very grateful to C. Keeton for a copy of 
his gravlens code and to the anonymous referee for a detailed veri- 
fication of our calculations which substantially improved the paper. 
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